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Abstract 



We have investigated analitycally the phase diagram of a generalized spher- 
ical version of the Blume-Emery-Griffiths model that includes ferromagnetic 
or antiferromagnetic spin interactions as well as quadrupole interactions in 
zero and nonzero magnetic field. We show that in three dimensions and 
zero magnetic field a regular paramagnetic-ferromagnetic (PM-FM) or a 
paramagnetic-antiferromagnetic (PM-AFM) phase transition occurs when- 
ever the magnetic spin interactions dominate over the quadrupole inter- 
actions. However, when spin and quadrupole interactions are important, 
there appears a reentrant FM-PM or AFM-PM phase transition at low 
temperatures, in addition to the regular PM-FM or PM-AFM phase tran- 
sitions. On the other hand, in a nonzero homogeneous external magnetic 
field H, we find no evidence of a transition to the state with spontaneous 
magnetization for FM interactions in three dimensions. Nonethelesss, for 
AFM interactions we do get a scenario similar to that described above for 
zero external magnetic field, except that the critical temperatures are now 
functions of H. We also find two critical field values, H c \, at which the 
reentrance phenomenon dissapears and H C 2 (H c i ~ 0.5-ff C 2), above which 
the PM-AFM transition temperature vanishes. 
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I. INTRODUCTION 



Continuous spins models play an important role as approximate statistical mechani- 
cal models for several physical systems. In particular, different versions of the classical 
spherical model (CSM) introduced by Berlin and Kac can be solved exactly for both, 
short [0 and long range |2||| interactions. The spherical model (SM) was originally in- 
troduced to study ferromagnetic systems in different spatial and spin dimensions fl]||. 
Unlike the mean-field approximation, in the case of spin dimensionality higher than one, 
the spherical model rightly yielded no ferromagnetic transition as had been observed ex- 
perimentally. Moreover, its equivalence to a spin system with infinity spin dimensions has 
been completly established by Stanley 0]. Because of this, different versions of the SM 
have been recently used to investigate the physics of systems which hamiltonian can be 
maped out onto that of the spherical model The basic idea is the compactification 

and replacing of too many restrictions by one spherical like condition. Thus, studying the 
phase transitions and the phase diagram of a spin ferromagnet or antiferromagnet using 
a particular version of the SM is of physical relevance. 

In this paper we investigate the phase behavior of a generalized spherical version of 
the Blume- Emery- Griffiths (BEG) model Q. This is obtained from the discrete three 
states spin one ferromagnet or antiferromagnet and, as in the original BEG model, we 
include magnetic spin interactions as well as quadrupole interactions. The BEG model 
was introduced to understand qualitatively the phase separation, and the A transition of 
3 He- 4 He mixtures. It has also been applied to study magnetic systems with competing 
interactions || or to describe the qualitative phase behavior of a microemulsion || . 

By introducing a transformation that maps out the original BEG model hamiltonian 
onto that of the SM we are able to investigate the conditions under which phase transitions 
may occur. Proceeding as in the investigation of the phase transitions of the CSM we 
study the properties of the saddle points of the integrand of the partition function. 

The main findings of this paper are the following: (1) There is no phase transition 
at all in d = 1, 2 spatial dimensions for both, the FM and AFM systems. This result 
agrees well with that obtained for the CSM. (2) For d = 3 dimensions and zero external 
magnetic field, we obtain the usual paramagnetic-ferromagnetic(antiferromagnetic) [PM- 
FM(AFM)] phase transition if the FM (AFM) interactions dominate over the quadrupole 
interactions. On the opposite side, when the quadrupole interactions dominate over the 
magnetic exchange couplings no phase transition occurs. This result is consistent with 
the limit of zero magnetic exchange interactions. However, in the intermediate regime, 
when both, quadrupole and exchange interactions are relevant, there appears at low tem- 
peratures, a novel and intriguing reentrant phase transition, in addition to the regular 
high temperature one. That is, the FM or AFM ordering takes place at intermediate 
temperatures, since at high and low temperatures a PM phase sets in. (3) For a FM 
system in a nonzero external uniform magnetic field and d = 3 spatial dimensions, the 
transition to a spontaneously magnetized state is destroyed, so no phase transition occurs. 
Nonetheless, for AFM exchange interactions in a uniform magnetic field we find a scenario 
similar to that described in (2) for H = 0, except that now the critical temperatures are 
now functions of H. More importantly, the domain of parameters where the reentrant 
phase transition happens narrows as the external magnetic field increases and, at the end, 
the transition completely disappears above a certain critical value H c \. Increasing further 



the external magnetic field yields a second critical field value H C 2 ~ 2H c i, above which 
the AFM phase is supressed. 

The layout of the remainder of this paper is the following: in section we introduce 
the BEG model hamiltonian and map it onto a generalized version of the spherical model. 
Then, we investigate the general conditions under which phase transitions may take place. 



In section |TTT] we consider the case of short range interactions and study the different 
scenarios that may occur in the system. A novel global spherical condition that is fully 
equivalent to the SM and valid even for low spin dimensionality is introduced in section 
|TV| . Finally, in section [V] we present a summary of our findings and the conclusions. 

II. THE MODEL 

Let us start with a BEG type model hamiltonian with spins = 1 located at each site 
of a d- dimensional rectangular lattice. Thus, each spin has three different states namely, 
= 0, ±1. It is natural to assume that the spins interact with each other through the 

pair potential V(R — R!) and with an effective external magnetic field H = MqHq, where 
Hq is a "true external field" and Mo is the magnetic moment produced by the states whith 
= ±1. According to these statements the Hamiltonian that defines the model is: 

w = ^£4 + £ v (£- s r s r' + # (i) 

R R,R! R 

where the second sum may run over all nearest neighbors or over all pairs of spins. The first 
term in 7i has been introduced in order to have the possibility of quadrupolar ordering, 
as suggested in reference [§J, in addition to the magnetic ordering associated with the 
second term. In absence of an external magnetic field the competition between these two 
types of ordering is measured by the parameter E, which can be interpreted as the energy 
difference between the states with = ±1 and = 0. With the aim at mapping this 
model hamiltonian onto that of the spherical model, we introduce the Ising-like variables 
d£ — ±1 and bg — ±1 that are related to the original spin variables by, 

Sr = \ {a & + b n ) . (2) 

Under this transformation the state Sg = has been double weighted since it can now 
be obtained from the states = —bg = 1 and = —b^ — — 1. To compensate 
this "degeneration" we introduce, in an ad hoc manner, the additional entropy term, 
Tis = —^T\ti2J2rS'^, in units in which the Boltzmann constant ks is equal to unity. In 
terms of these Ising-like new variables and including the entropy term the hamiltonian TC 
becomes: 

« = \( E ~ Tln2 ) £ (a M + b n ) 2 + ±Y J V(R-R') {a A + b A ) (a M , + b &l ) 

R R,R' 

+ f E( fl a + 6 s)- ( 3 ) 

R 

To be able to study analytically the properties of the model we have to make contact with 
the classical spherical model So, we replace the discrete Ising- variables by continuous 



ones, that is, we let the "new" spin variables to vary in the interval, — oo < ag, < +00, 
obeying the spherical conditions: 

£4 = A/", and (4) 

R R 

where M is the total number of spins in the system. Hence, the hamiltonian, Eqn. (|3]) 
together with the last two restrictions can be considered as the "generalized spherical ver- 
sion" of the BEG hamiltonian, Eqn. (|l]). The reason of this name is that we have included 
the additional entropy term that depends explicitly on temperature in the quadrupole in- 
teraction part of the Hamiltonian. This term does not appear in the original BEG model 
and will play a key role in the results we describe in what follows. 

Taking into account these spherical conditions on the "new" Ising-like variables, one 
can write the partition function as follows: 

2 = [ n £L da A SZ ^ j ex p 5 - e 4) 5 - e 4) • ( 5 ) 

It is important to point out that the transformation from discrete Ising spin variables 
to the continuous spin ones is not quite correct since there are "wrong" configurations in 
which there are spin variables that may take values different from ±1. Nonetheless, such 
"slight incorrect transformation" allows one to solve analitycally two classical problems: 
the 3D Ising Model [0,0, and more importantly, the 3D Heisenberg model |7|, in addition 
to other more recent contemporary physical problems p|-pJ|. In section IV we address 
this shortcoming and propose an absolutely correct global spherical condition that is even 
valid for systems of low spin dimensionality, in particular for = 1. For the time being 
we proceed as in the CSM, that is, we use the integral representation of 5-function to 
rewrite the partition function as: 

Z = /|0exp(AT( Sl + S2 ))z( Sl , S2 ), (6) 



where 



2= {HI dag I (7) 



if 



—00 j —00 



and 



H = Ts 1 J2a% + Ts 2 J2b% + H, 

if 



where Tt is given in Eqn. (|1|), with the replaced by |(^ + ^)- in the framework of the 
CSM the effective hamiltonian 7i is a non-diagonal symmetric bilinear form in the new 
Ising-like variables. To be able to integrate Eqn. (0) one has to diagonalize 7i by means 
of a Fourier transformation of the spin variables. In doing so one gets, 



where 



Si + So t? ^ In 

IN ^ 



SlS 2 + (si + s 2 



and 



T 



+ 



H 2 



81 + s 2 



lQrp2 o(0) 
T 



Si + S 2 ) + SiS 2 



(10) 



E -T\n2 + J2V (r) exp(-iqR) 



R 



where the summation in Eqn. fllCf ) is carried out over the first Brillouin zone (BZ). The 
next step would be the evaluation of the partition function, however, instead of doing 
it we proceed as in the study of the CSM. That is, we find the saddle points of Z 
and study its temperature behavior since it contains the complete information of the 
possible phase transitions of the system. Nonetheless, it is important to note that the 
evaluation of the partition function by the steepest descendent method becomes exact in 
the thermodynamic limit M —>■ 00, as in the CSM. The saddle points of the integrand of 
Z are determined by the following equations: 



dsi 



0=1- 



H 2 



1 



2( Sl + s 2 ) 27Vsi + s 2 



E 



1 



Sl S 2 + ^ (Si + s 2 ) 



16T 2 



(si + 8 2 ) + SiS 2 



2 ' 



(12) 



^l = o = i I ^^L_y I 

ds 2 2( Sl + s 2 ) 2Af Sl + s 2 ^ SlS2 + ^( Sl + s 2 ) 

H 2 s 2 

(13) 



16T2 ( t 1 («i + *>) + s,s 2 

These equations are symmetric in the variables s\ and s 2 . Substracting Eqn. flllf ) from 
Eqn. ( |T^ ) and taking into account that 7(^1, s 2 ) has to be an analytic function at the 
saddle point, one gets the condition si = s 2 = s. Substituting this in either, Eqns. ([12]) or 
(|13|), we obtain the equation that yields the saddle points of the integrand of the partition 
function, namely 

*>4?7tV%^ + ; =4 ' (14) 

To investigate the existence of phase transitions in the model, we should recall that, 
as in the CSM, the existence of a saddle point is associated with the setting in of a 
disordered phase, PM phase in this case, while the absence of a saddle point is related 
to the realization of an ordered phase, FM or AFM. The critical temperature is obtained 
by studying the evolution of the saddle point as a function of temperature. To be more 



precise, the phase transition occurs just when the saddle point dissapears. As we will 
see below there is a such transition in some region of the parameters space (E,J,H), 
nonetheless, there is another region in the parameters space where a saddle point does 
exist for any temperature, indicating the absence of a phase transition in such region. 
For instance, one can easily show that for a nonzero external magnetic field and if the 
ferromagnetic interactions between spins satisfy the condition mm^a(q) = a(0), then 
a saddle point exists in the whole parameter space and therefore no phase transition 
takes place. However, for antiferromagnetic interactions one gets min ^a(q) < a(0) and 
a more careful and detailed analysis should be carried out to obtain the corresponding 
phase diagram. (See subsection [111 B| ). In the following section we apply our general model 
hamiltonian to study a system with short range interactions. 



III. PHASE TRANSITIONS IN A FERROMAGNETIC AND AN 
ANTIFERROMAGNETIC SYSTEMS WITH SHORT-RANGE 

INTERACTIONS. 



A. Zero external magnetic field. 

Let us consider a magnetic system with short range interactions, that is, V(R — R') 
—J for nearest neighbors R, R', and zero otherwise. In this case Eqn. (|Il]) becomes: 



1 



d 



a(cfi = -\E-Tlri2-2jJ2cosq j J, (15) 

where d is the spatial dimensionality, q = (ai<7i, 0,20.2 ■, ) is a dimensionless vector in the 

reciprocal space and a, is the lattice constant along the j direction. Setting H = in 
Eqn. (|T4]), the function if) (s) = ipo (s) can be rewritten as 

* = = ; + j<k?L d « 5 ( s+ |- |h !)- ( Etl cos „) • (16) 

here the integration is carried out in the first Brillouin zone, that is, in the range — 7r < 
qj < 7r for each j. As we know and according to the definition of the interaction potential, 
J > 0, for ferromagnetic interactions while J < 0, for antiferromagnetic interactions. 
Notice that in the framework of our model the phase diagram in the E versus J plane 
for the antiferromagnet (J < 0) is obtained from the reflection about the line J = 
of the phase diagram of the ferromagnet. Indeed, by introducing a shift in the wave 
vectors in g-space, that is, by replacing q by q* = (q\ — 7r, g 2 _ tt, •••) one can show that 
JJ2 cos qj = —JJ^cosq*. Therefore, Eqn. (|T6|) for the antiferromagnetic system results in 
that corresponding to a ferromagnet. The only difference being the critical value q = q m 
at which a(q) attains its minimum value. For a ferromagnet this minimum is reached at 
q m = while for an antiferromagnet the minimum occurs at q m = (— ir, —it, ...). In term of 
the orientations of spins this is equivalent to say that nearest-neighbors spins try to align 
parallel for ferromagnetic interactions or anti-parallel for antiferromagnetic interactions. 
Thus, all particularities of the phase transition of the antiferromagnet are the same as 
those for the ferromagnet. Bearing this in mind we carry out a detailed analysis of the 



phase diagram for a ferromagnet and extend directly the results for an antiferromagnet 
by carrying out a reflection about J = 0. 

The next step is to examine the properties of i/jo(s) for a ferromagnet since they 
determine the phase transitions in the model under investigation. It is easy to see that 
V'o(s) decreases monotonically and is analytic in the interval Sq(T) < s < +00, where so 
is defined as, 

s = max ^0, - In 2 - (E - 2d J) /2T\ . (17) 

The behavior of the function "0o( s ) is such that it can take values from zero when s — ► 00 
up to a nonzero value t/>o( s o)- Its behavior is shown in Fig. [TJ, where we have indicated 
the positions of Sq (points B and A) at two temperatures T\ and T 2 that are below and 
above the critical temperature T c , respectively. The point R represents the saddle point 
which existence is related to the set in of the disordered phase. Thus, at T = T 2 (point P 



and dot dashed line), ^o( s o) > 4, and Eqn. ( [Ty) has a root, indicative of the existence of 



a saddle point and a paramagnetic phase. In the opposite case, at T — T\ (point Q and 
solid line), Vo(so) < 4, and Eqn. ( |l(| ) has no root yielding no phase transition to the 
paramagnetic phase. Hence it is necesary to analyze the behavior of tpo(so) as a function 
of temperature in the parameters space (E, J) to be able to construct the phase diagram. 
First of all, we have to study the general conditions for which the inequality ipo( s o) < 4, 
responsible for the ferromagnetic phase, is valid. Thus, in the region of parameters (E, 
J,T) where s = 0, the function ipo( s o) — * 00 (^o( s o) > 4) and a paramagnetic phase 
exists. Nonetheless, when so > it is necessary to analyse the convergence of the integral 
for s = so in Eqn. (|16|). In doing so we substitute s = so and expand cosg^ ~ 1 — g|/2, 
yielding, 

1 f ^dq (w) 



*z. 5( So + J, _ I l n 2 ) -(£? = i cos?,-) J b.z. 

where we have omited an irrelevant constant numerical factor. The last integral converges 
if the spatial dimension d is such that, d — 3 > —1, i.e. for spatial dimensionality d > 2. 
This means that for d < 2, tpo(s — > s ) — > 00 yielding a saddle point for any values of 
E, J and T, and no phase transition takes place for lower dimensional systems. A more 
detailed analysis has to be carried out for three dimensions because the integral fli~8|) does 
converge and can be evaluated as Jl], 

If dq 

' ' 0.50541 = L r . 



(2tt) 3 7 3-E^iCosg, 
Introducing this result in Eqn. fllBD one obtains, for Sq(T) > 



T 2T 

<M«<r>,T) = 1 i. + Tru2 _ E+tJ . (19) 



In this case we find (see Fig. 0) that there are regions in the parameter space (E, J, T) 
where ^o( s o) < 4, i.e. the saddle point disappears. This happens for temperatures 
below a certain critical temperature, T c , indicating the transition to the ferromagnetic 



(antiferromagnetic) state. A more detailed analysis of ■?/>o( s o(^)> Tj shows that there are, 
actually, two ranges of parameters to be considered: (1) E — 6J < in which case s (T) 
goes from +00 at T — up to |fn2 at T — 00, and (2) E — 6 J > for which s (T) goes 
from zero at T = up to | In 2 at T = 00. 

Let us first consider the behavior of i(jq(sq(T),T) as a function of temperature when 
E — 6 J < 0. In this domain of parameters Sq > for any T and Eqn. ( |T9D is always 
valid. At very low temperatures ip (s (T),T) < 4 and at T = 0, ipo{ s o) g° es down to zero, 
implying that there is no saddle point, thus the system settles down in a ferromagnetic 
state. At higher temperatures, the function ip (so{T),T) increases almost linearly with 
temperature and therefore ^o( s o) crosses over the line ipo( s o) — 4, at the critical tem- 
perature and the saddle point shows up, hence the system goes from the ferromagnetic 
(antiferromagnetic) phase to the paramagnetic one. Thus, this region of parameters yields 
one critical temperature that corresponds to the usual PM-FM (or PM-AFM) phase tran- 
sition. On the other hand, in the range of parameters where E — 6 J > the behavior of 
ijjo(so(T), T) is different. For T < (E — 6J)/ ln2, the function V'o(so) diverges since so = 0. 
For s (T) > 0, the function ip (s (T)) decreases asymptotically and quite rapidly from 
00 at T = (E — 6 J)/ In 2 until it reaches a minimum at the temperature T* = T*(E, J), 
given by 



1 cr 



and afterwards, it begins to increase. This behavior is shown in detail in Fig. 0. It is 
obvious that as soon as -?/>o( s o(^*)) T*) > 4 then i/; (s (T),T) > 4 for any temperature 
T greater than T* and the paramagnetic phase is only present, therefore, no phase tran- 
sition occurs. This is what happens when j > 6 + x c for any value of T, where, the 
value of x c is given by x c = (l// cr )(2v / hi2 - V2) 2 w 0.124 . The corresponding curves 
ij)o(so(T),T) that fullfil this condition are shown in Fig. |2|. However, if 6 < j < 6 + x c 
then ip (s (T*),T*) < 4, and ^ ( s o(^); T) intersects the line ip = 4 in two points, yielding 
two critical temperatures: T c + , the branch where the PM-FM (PM-AFM) phase transition 
occurs, and T~, the branch at which the system goes back to the paramagnetic phase; so 
a reentrant phase transition takes place. These critical temperatures are given by 

rr± _ (E-QJ)I cv + AJ\n2-2J±V 

1 c — or i„o ' y ZL ) 



with 



V 



2I CT In 2 



((E — 6J)/ cr + 4 Jin 2 — 2 J) 2 — 16(£ - 6 J)/ cr J In 2 (22) 



1/2 



From these results we can state that in the temperature range T~ < T < T c + a magnetic 
ordered state takes place, whereas outside of this interval, i.e. at lower T < T~ and 
higher T > T c + temperatures, a paramagnetic state settles in. Translating these results 
in the parameter domain we can say that, the region 6 < -j < 6 + x c where the reentrant 
phase transition occurs, separates out the region E — 6J < where the usual PM-FM 
(PM-AFM) phase transition happens, from the region 6 + x c < y where there is no phase 
transition at all. 



To make more intuitive and explicit the above results we have plotted the two dime- 
sional phase diagram ^j- versus -j in figure The solid line represents the critical 
temperature T c + that corresponds to the usual high temperature single phase transition 
while the dashed line represents the novel and intriguing low temperature reentrant phase 
transition T~ . These lines have also been plotted in a larger scale in the inset, where the 
two branches are clearly seen. The two critical temperatures become gradually closer to 
each other with increasing E/J until the ordered phase dissapears completely. 



B. Nonzero uniform magnetic field. 

In this subsection we investigate how the phase diagram of the generalized BEG model 
changes in the presence of an homogeneous magnetic field in three spatial dimensions. 
In the case of ferromagnetic interactions the application of an external magnetic field 
terminates the transition to the state with a spontaneous magnetization since a saddle 
point appears for any values of E, J and T, as it just happens for the CSM Jl|. That 
is, the external magnetic field is a symmetry breaking field that produces a magnetized 
state that always yields a saddle point in the integrand of the partition function for any 
values of the parameters. Because of this, an additional analysis should be done to study 
the thermodynamics of the present model, for instance, one can study the behavior of 
the magnetic susceptibility, as was done in |IJ] for the CSM. For the time being, it is 
the purpose of the present paper to study the phase diagram and leave the study of the 
thermodynamics for a forthcoming paper Let us now focus our attention on the 

case of antiferromagnetic interactions. Basically, we will study the circumstances under 
which the saddle point dissapears in the presence of an homogeneous magnetic field, i.e. 
when the antiferromagnetic ordering sets in. Using Eqn. (|i4|) and the results of section 



I [I A| it is easy to show that ip(s) is an analytic function that decreases monotonically in 
the interval s < s < oo, where s is given by Eqn. ([T7|). To be able to use the results 
obtained for an antiferromagnet in a zero magnetic field, in what follows we will denote 
the modulus of the AFM exchange interaction as J > since the minus sign has already 
been included as a shift in the reciprocal space vectors, that is, we have replaced q by q*. 
Henceforth, in analogy with the case of zero magnetic field, whenever ip(so(T),T, H) > 4, 
a saddle point does exist and it yields a paramagnetic phase, whereas in the opposite case, 
i/j(sq(T),T, H) < 4, the system settles in an antiferromagnetic state, where the function 
ip(s (T),T, H) is given by 

i, { s (T),T,H) = M So(T),T) + (£jf = + TlB2 2T E + 6J + (23) 

To obtain this equation we have made the following substitution in Eqn. ([14]), 

tf 2 /(4T 2 ( So + 2a(q = 0)/T) 2 ) = {^f . 

As in the previous subsection, the equation i/jq(sq(T c ),T c , H) = 4 determines the temper- 
ature at which the phase transition occurs. From this very last condition one can rewrite 
Eqn. (|23|) as , 

^ ( So (T c ),T c )=4-( I | 7 ) 2 , (24) 



where ^q(sq(T),T) is the function that we have carefully studied for the ferromagnetic and 
antiferromagnetic systems in abscence of an external field. Thus, one obtains T C (H) by 
looking at the intercept of the curves ipoi s o) shown in Fig. |2|, with the lines 4 — (if/12 J) 2 . 
The resulting plots are shown in Fig. (| for the following values of the reduced magnetic 
field, f : (a) 2, (b) 8, (c) 16 and (d) 20. One immediately sees that in the parameter region 
E — 6 J < 0, there is only one root of the equation ip(so(T c ), T c , H) = 4, what signals the 
usual transition to the antiferromagnetic phase. The critical temperature of this transi- 
tion, T C (H), that is now a function of the external magnetic field, diminishes as the field 
intensity increases and becomes equal to zero at the critical field value H c2 = 24J. Thus, 
the phase transition dissapears for fields such that H > H c2 , and the paramagnetic state 
settles in at all temperatures, as one naturally would expect, since the antiferromagnetic 
ordering is broken by higher external field intensities. We have represented this critical 
field value by H C 2 because of, as we will show below, there is another critical field value 
H c i fa 0.5H C 2 at which the reentrant phase transition dissapears. 

On the other hand, in the parameter region where E — 6 J > there are two possibili- 
ties: (i) the reentrant AFM-PM phase transition or (ii) the absence of a phase transition 
at all. The domain of parameters where the reentrant transition takes place is restricted 
by the inequalities, E - 6 J > and vjj (s (T*),T*) < 4 - (H/12J) 2 , where T* is defined 
in Eqn. (p0|). These two relations can be reexpressed in the form: 6 < Ej J < 6 + x c (H), 
with x c (H) given by 



2Vhi2,/l - (^) 2 - y/2 



x c = ± * '-. (25) 

-*cr 

The two critical temperature branches of the phase transitions are the roots of the equation 
ip(so(T c ),T c , H) = 4, and are given by 



(E - 6J)/ cr -2J + JU- (H/12 J) 2 In 2 ± V(H) 

= kh2 • < 26 > 



where 



V = \J ((6J - E)I CI + 2J - J (4 - (^jf) ln2) + 4/ cr J(6J - E) (4 - (^) 2 ) h, 2. 



(27) 



Evidently, the reentrance phenomenon is suppressed by the external field when it exceeds 
the critical value H c \ defined as the root of the equation x c (H c i) = 0. This leads to the 
critical field value, 



'4-2/ln2 Ho 
H cl = H c2 ^ J- = 0.52788# c2 w -f. (28) 

These results allow us to construct the phase diagram (j versus 7, y) which, for the sake 
of clarity, is shown in two parts, Figs. |5| and |^. The former shows the region in three 
dimensional space where there is a single transition, while the latter corresponds to the 
region where the reentrance phenomenon takes place. That is, the whole half-plane ( j, 7) 



splits into three regions, the first one corresponds to the "single transition" characterized 
by the PM-AFM phase transition, while the "reentrant AFM-PM phase transition" occurs 
in the second region. The third region corresponds to the domain where there is no phase 
transition at any temperature and any values of the coupling parameters E and J. 



IV. EXACT GLOBAL SPHERICAL CONDITION. 

One of the main shortcomings of SM approach, that we have used here, Eqn. (f|), is 
that it fails when spatial and spin dimensionality are low. For instance, it is known that 
in two dimensions the CSM approximation for an Ising ferromagnet does not describe 
the correct transition scenario. Nonetheless, for spin dimensions greater than one, the 
CSM approximation works quite well. In fact, it is exact in the limit of infinite spin 
dimensionality |7]]. The main drawback of the classical spherical condition: Y^ni a \~ 1) = 
is that it holds not only for "right configurations" , = ±1, but it also does for the "wrong 
configurations" where some of the spin variables are < 1 while others are ag> 1. 

So, in trying to fix up this shortcoming we replace the usual spherical condition by 
the "general global" one, namely 

E(4 - !) 2 = 0- (29) 

R 

This sum is equal to zero only for the "right configurations" = ±1. Therefore, the 
model with continuously varying spin variables plus the "global spherical condition", 



Eqn. (29), are equivalent to the original spin one BEG hamiltonian. Note that we may 
also use the general global condition J2jif( a \ — 1) = where f(x) > for x ^ and 
/(0) = 0, however, for definiteness and simplicity one can use Eqn. (p9|). This condition 



can be expressed in a form slightly closer to the usual spherical condition as follows: 

-^ = E4 + E(4- 1 )(4-2)- (3°) 



R R 



With this "new global condition" we can write a new effective hamiltonian which in 
turn defines the partition function Z(si,s 2 ) (see, Eqns. (JJ) and (§)): 

H( Sl , s 2 ) =H + T Sl 5>| - l)(a| - 2) + Ts 2 E(&| - 1)(&| - 2), (31) 

R R 

where Tio is the usual spherical effective hamiltonian defined by Eqn. (§). So one can 
rewrite the new effective hamiltonian as the usual spherical hamiltonian plus an additional 
interaction term U(si, s 2 ), that is: H(si, s 2 ) = Ho + U(si, s 2 ). Note that this new effective 
hamiltonian is much more complicated than the one we have studied in this paper and, 
unfortunately, the interaction [/(si,s 2 ) cannot be considered perturbatively because it 
lacks of a small parameter. However, our propossal allows us, in principle, to improve the 
classical spherical model for a system with lower spin and space dimensionalities. The 
problem of incorporating the U(s) interaction into the complete solution is a difficult task 
and it is the subject of current research. 



V. SUMMARY AND CONCLUSIONS 



In this paper we have introduced and studied the phase diagrams of a generalized 
spherical version of the BEG model considering ferromagnetic as well as antiferromag- 
netic interactions. The model hamiltonian involves terms representative of quadrupole 
interactions and magnetic spin interactions in zero and nonzero external magnetic field. 
We have shown explicitly that in the short range interactions limit and zero external 
magnetic field the model presents no phase transition in one and two spatial dimensions. 
However, in three dimensions it undergoes a regular PM-FM(AFM) transition at nonzero 
temperature. This happens in the case when the magnetic interactions dominate over 
the quadrupole ones. Nonetheless, in the range of parameters where quadrupole and 
magnetic interactions are relevant we obtain a low temperature reentrant FM(AFM)-PM 
phase transition, in addition to the conventional PM-FM(AFM) phase transition at higher 
temperatures. These phase transitions are absent in the strong quadrupole interaction 
regime. On the other hand, when there is an external magnetic field there is no phase 
transition for the ferromagnet in three dimensions, since there appears a permanent mag- 
netized state that terminates the PM-FM transition. This result appears to be similar 
to the one obtained in the CSM 0, and it is a motivation to further investigate the 
thermodynamics of the present model However, in the case of antiferromagnetic 



interactions and nonzero magnetic field we also get the reentance phenomenon whenever 
the magnetic spin and quadrupole interactions compete between them and always that 
H < H c \ 12J. If H c i < H the reentrant phase transition dissapears. We also find a 
second critical value of the external magnetic field above which the PM-AFM phase 
transition dissapears. Finally, we have introduced a "novel global spherical condition" 
that makes our approach valid even for low spin dimensionality. However, this led us to 
a more complicated hamiltonian which phase diagram is under current investigation. 
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FIG. 1. Behavior of ip as function of the parameter s at two temperatures T\ and T2 below 
and above the critical temperature, respectively. Points A and B have coordinates ^so(2~2), O^j 

and ^so(^i)>o)j while the coordinates of points P and Q are ^0(^2),^ and (so(Ti) , ip*j , 
respectively. Point R indicates the saddle point that signals the phase transition. Note that ip(s) 
decreases monotonically with s, that is, it goes from Y>(so) and goes down to zero asymptotically 
as s — > 00. 

FIG. 2. Behavior of Y>o(so) as function of j and different values of -j. Note that for 
< -j < 6 there is only one intercept of these curves with the line Vo = 4, as an indication 
of a single transition. However, for 6 < -j < 6 + x c we see a lower temperature intercept 
corresponding to the reentrant phase transition. 



FIG. 3. Phase diagram for the ferromagnet where the reentrant transition appears at -j > 6. 
The letters P and F indicate the paramagnetic and ferromagnetic states, respectively. In the 
inset we show the part of the phase diagram corresponding to the reentrance behavior in the 
parameter region 6 < -j < 6 + x c . 



FIG. 4. Behavior of ipo(so) versus -j for different values of the parameter -j (same values 
as in Fig. |2|) and different values of the external magnetic field, (a) -j = 2 (b) -j = 8, (c) 
f = 16, and (d) f = 20. Notice that for E/J < 6, where a single transition occurs, the critical 
temperatures T C (H), decrease as the magnetic field increases. On the other hand, for curves 
ip(s), with 6 < -j < 6 + x c (H), there are two transitions, the high temperature one and the one 
corresponding to the reentrance to the paramagnetic state, panels (a) and (b), while in panels 
(c) and (d) there is no phase transition at all, because of the fact that -j > £f. 



FIG. 5. Three dimensional phase diagram -y versus -j, y in the range of parameters where 
there is a single transition. 



FIG. 6. Three dimensional phase diagram versus -j, -j restricted to the range of param- 
eters where the reentrant transition takes place. 
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Fig. 5 (S.E. Savel'ev and G. Ramirez-Santiago) 



